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Abstract. Respiration-correlated CBCT, commonly called 4DCBCT, provide 
respiratory phase-resolved CBCT images. In many clinical applications, it is more 
preferable to reconstruct true 4DCBCT with the 4th dimension being time, i.e., each 
CBCT image is reconstructed based on the corresponding instantaneous projection. 
We propose in this work a novel algorithm for the reconstruction of this truly time- 
resolved CBCT, called cine-CBCT, by effectively utilizing the underlying temporal 
coherence, such as periodicity or repetition, in those cine-CBCT images. Assuming 
each column of the matrix U represents a CBCT image to be reconstructed and the 
total number of columns is the same as the number of projections, the central idea of 
our algorithm is that the rank of U is much smaller than the number of projections 
and we can use a matrix factorization form U = LR for U. The number of columns 
for the matrix L constraints the rank of U and hence implicitly imposing a temporal 
coherence condition among all the images in cine-CBCT. The desired image properties 
in L and the periodicity of the breathing pattern are achieved by penalizing the sparsity 
of the tight wavelet frame transform of L and that of the Fourier transform of R, 
respectively. A split Brcgman method is used to solve the problem. In this paper 
we focus on presenting this new algorithm and showing the proof of principle using 
simulation studies on an NCAT phantom. 
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1. Introduction 

When cone beam computed tomography (CBCT) is applied to thorax or upper abdomen 
regions, motion-induced artifacts, such as blurring or distortion, become a serious 
problem, because different x-ray projections correspond to different volumetric CBCT 
images due to patient respiratory motion. To overcome this problem, four-dimensional 
CBCT (4DCBCT)(Sonke et al, 2005; Kriminski et al, 2005; Li et al., 2006) has been 
developed. In such a modality, all x-ray projections are first retrospectively grouped 
into different respiratory phase bins according to a breathing signal tagged on every 
projection image. A set of CBCT images are then reconstructed, each at a breathing 
phase, under the assumption that the projections placed into each bin correspond to 
the same or similar CBCT image. The number of phase bins is usually empirically 
chosen by a user based on the consideration of balancing temporal resolution and image 
quality. On one hand, a high temporal resolution requires a large number of phase 
bins, which leads to insufficient number of projections available to each phase and hence 
degraded quality of reconstructed CBCT images. On the other, a small number of phase 
bins results in low temporal resolution, as well as a relatively large bin width and the 
associated residual motion artifacts in the reconstructed images. To maintain clinically 
acceptable temporal resolution and image quality, 4DCBCT acquisition protocols such 
as slow gantry rotations and multiple gantry rotations have been invented to increase 
the total number of projections(Li et al., 2006; Lu et al., 2007; Li and Xing, 2007). 

Conventional 4DCBCT reconstruction approach reconstructs CBCT images at 
different phases independent of each other (Sonke et al., 2005). This straightforward 
method neglects the temporal correlations of CBCT images at different phases. 
In contrast, a collaborative reconstruction scheme has been recently proposed(Jia 
et al., 20f0; Gao et al., 2011; Tian et al., 2011; Qi and Chen, 2011; Jia et al, 2012). 
Among them, Jia et al. (2010) and Tian et al. (2011) proposed a temporal non-local 
means method to constraint that the reconstructed CBCT images at neighboring phases 
must contain repetitive anatomical features. The locations of these features are allowed 
to vary among phases. Qi and Chen (2011) first reconstruct an average CBCT image 
using all projections at all phases and impose a similarity constraint between this 
average image and the CBCT at each phase. These methods explicitly enforces the 
similarity among reconstructed CBCT images. Gao et al. (2011) utilizes robust PCA 
techniques for the 4DCBCT reconstruction problem. This method restores a matrix 
whose columns are the 4DCBCT images at all phases under the assumption that, after 
a proper transformation, this matrix can be decomposed as the sum of a low-rank matrix 
corresponding to an almost static background and a sparse matrix corresponding to a 
moving foreground. The low-rank constraint achieved by minimizing the nuclear norm 
(Cai et al., 2010) of the associated matrix implicitly imposes the inter-phase similarity 
of the background. However, in the context of anatomical motion in thorax or upper 
abdomen, there is no clear distinction between foreground and background. Besides, 
as pointed out in (Chandrasekaran et al., 2010), the nuclear norm minimization always 
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finds a matrix factorization with mutually orthogonal factors. In this sense, it may have 
some unfavourable bias on the resulting low-rank matrices. 

Although commonly called 4DCBCT, respiration-correlated CBCT (Ford et al., 
2003) is a more accurate name for this imaging technique, since its 4th dimension is 
actually the respiratory phase rather than time. Clinically, it is more preferable to 
reconstruct true 4DCBCT with the 4th dimension being time, i.e., one CBCT image 
is reconstructed based on the corresponding instantaneous projection, if possible. This 
truly time-resolved CBCT, called cine-CBCT in this work to avoid confusion, encounters 
an apparent technical challenge of reconstructing a CBCT image based on only one 
projection. Previously we have tried to overcome this challenge by utilizing prior 
4DCT images of the same patient (Li et al., 2010; Li et al., 2011). Particularly, we 
have built a lung motion model for the patient by performing a principle component 
analysis(PCA) on the motion vector fields obtained from deformable registration on 
prior 4DCT images. We discovered that only a few principle components are sufficient 
to represent the lung motion to a satisfactory degree of accuracy. The reconstruction of 
a volumetric image based on one acquired x-ray projection is then achieved by finding 
those PCA coefficients, so that the projection of the corresponding CBCT matches the 
acquired one. However, this method heavily depends on how similar the current CBCT 
images are to the prior 4DCT images and how accurately the PCA parameters can be 
determined using one projection. 

In this paper, we will show the principle of proof of an innovative algorithm for cine- 
CBCT reconstruction without using any prior images. The basic idea is to maximally 
exploit the temporal correlations, both local and nonlocal, of the patient anatomy at 
various time points, namely smooth variation of CBCT images at successive time points 
and the periodicity of all images during a scan, using recent mathematical developments. 
Specifically, a factorization of the matrix consisting of all the CBCT images into a 
product of two matrices is used to explicitly enforce a low-rank condition of the matrix 
that is implied by the local or nonlocal temporal image similarity. Moreover, sparsity 
conditions are enforced under appropriate transforms, i.e. wavelet transform in space 
for one matrix and Fourier transform in time for the other matrix, to satisfy desired 
properties for the sequence of CBCT images in both space and time. Simulation studies 
have demonstrated promising results for our method for this challenging problem. 

2. Methods and Materials 

2.1. Model 

In this proof-of-principle study, we consider the reconstruction of a 2D slice of the CBCT 
to illustrate the principles of the new algorithm. The basic ideas, however, can be easily 
generalized into 3D contexts. Let the unknown 2D image be a function u defined in R 2 . 
When the X-ray source is placed with angle 9, the projection measured by the imager 



Cine-CBCT reconstruction 



5 



at location z is 

rt, 

fe(z) = V e u{z) := / u(x e + sr z )ds (1) 
Jo 

where Xg G IR 2 is the coordinate of the X-ray source, and r z G M 2 and £ 2 are respectively 
the direction and the length of the line connecting the X-ray source and the location 
z on the imager. The operator V in ([1]) is also known as Radon transform. If fg(z) 
is sampled with respect to z, the resulting projection data can be written as a vector 
fe G R A/ that obeys 

Pen = f e . (2) 

Here Pg G ]R MxAf and u G M. N are the discretization of Vg and tt in (JTJ respectively. 
Suppose that we have T projections where the X-ray source is placed with angels 
6*i, 6*2, • • • , 6t respectively. Let Pi and fi stand for Pg i and fg i respectively. Then, 
by putting (j2J) with different angles together, the CT projection can be written into a 
system of linear equations 

Pu = /, (3) 

where P = [P x ; P 2 , P T ] G R MT * N and / = f 2 ; f T ] G R MT . In other 

words, the CT image reconstruction problem is to recover image u from its partial 
Radon transform. 

In the context of cine-CBCT reconstruction, instead of only one unknown image 
u as in the CBCT problem, there is a set of unknown images, each associated with 
a projection, denoted by {ui, u 2 , ■ ■ ■ , Ut} and the projection condition in (J3J) may be 
modified to 

P lUl = f u z = l,2,...,T. (4) 

If we write all the unknown images in a matrix form U = [u±, u 2 , ■ ■ ■ , Ut\, (SJ) can be 
rewritten into a compact matrix equation 

VU = F (5) 

where PU = [P lUl , . . . , P T u T ] and F = [fi, . . . , f T ]. 

Because of the temporal coherence among the unknown images, it is expected that 
the columns of the matrix U formulated as such, namely the CBCT images, are almost 
linear dependent, and hence the rank of the resulting matrix U is very low, much smaller 
than the number of probections. With the low rank assumption, the number of intrinsic 
unknowns of (jSJ) are possibly less than the number of measurements, and therefore it 
is possible to reconstruct T CBCT images in U from T projections. To incorporate 
the underlying low-rank assumption into the cine-CBCT reconstruction process, we 
would like to explore matrix factorization based low-rank models in this study. More 
specifically, we enforce a decomposition form of the unknown matrix U of images as 
U = LR where L G R Arxir and R G M^ xT for a small integer K. From basic linear 
algebra knowledge, any matrix U of rank K can be represented in U = LR and 
conversely, the rank of U is at most K given the factorization form. In our algorithm, 



Cine-CBCT reconstruction 



6 



the value of K is specified by the user and our method regarding the selection of K will 
be presented later. We would like to point out that there is another popular approach 
for imposing the low rank condition that penalizes the nuclear norm of the matrix U, 
namely the sum of the singular values. Compare these two methods, the latter always 
finds orthogonal basis because of the involved singular value decomposition process. 
In contrast, our factorization method does not require this implicit orthogonality and 
hence attains the advantage of avoiding unfavourable bias on the resulting low-rank 
matrices (Chandrasekaran et al., 2010). In some sense, the approach of penalizing 
the nuclear norm of U is like Principle Component Analysis and our approach is like 
Independent Component Analysis. As a consequence, each of our basis, i.e., each column 
of L, represents more meaningful images and hence we can enforce some desired image 
properties on L, such as sparsity under wavelet transform. 

This matrix factorization approach also allows for the additional regularizations on 
L and R, so that U = LR carries desirable physical properties. First, the columns 
of the matrix L, corresponding to images, can be interpreted as a basis that efficiently 
represents all the columns of U. Notice that images usually have sparse coefficients 
under suitable transforms such as wavelet tight frames (Ron and Shen, 1997). Let T> be 
such a transform. Since the columns of L are images, we want T>L to be sparse. It is 
well-known that ^-norm minimization leads to sparse solutions. Therefore, the ^-norm 
||X>.L||i is penalized. Second, the rows of R are the coefficients of U under the basis L. 
Since the respiratory motion is near periodic, each row of R can be modeled as a period 
signal in time. In other words, the Fourier transform of each row of R contains very 
few nonzero entries. Since the fx-norm minimization helps get sparse solutions, we can 
use ||J-\R||i as a regularization term for R, where J 7 is the Fourier transform applied to 
each row. 

Altogether, we propose to reconstruct the cine-CBCT by solving an optimization 
problem 

min||X>L||i + A||J 7 fl||i, s.t. V{LR) = F, (6) 

where A is a parameter that balances the sparsity of T>L and JFR. In practice, (jSJ) 
is barely satisfied, because there is always unavoidable noise in the measurements 
F. Moreover, in reality the matrix U is only approximately low-rank and explicitly 
enforcing a low-rank representation leads to error. For these considerations, we solve 

min||X>L||i + A||J 7 i2||i, s.t. \\V{LR) - F\\ 2 F < a 2 , (7) 

where || • \\p is the Frobenius norm, namely || A\\p = yj^ij A 2 - for a matrix A. In (jTJ), 
A is again a parameter balancing the sparsity of T>L and FR, and a is a parameter to 
control to what extent the violation of (jSJ) is allowed. 

2.2. Algorithm 

Let us first consider the algorithm for solving Equation (Q. We use a split Bregman 
method (Goldstein and Osher, 2009; Cai et al., 20096) (also known as augmented 
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Lagrangian method) to solve this problem. In particular, by introducing some auxiliary 
variables C and D, the problem in is equivalent to 

min ||C||i + AllDlli, s.t. V(LR) = F, C = VL, D = TR. (8) 

C,D,L,R 

The augmented Lagrangian is 

E(C, D, L, R, Z, Zi, Z 2 ) = ||C||i + All-Dili + (Z, V(LR) - F) + |||7>(M) - F|£ 

+(Z 1; C - Pi) + ^||C - PL||| + (Z 2 , D-FR) + ^\\D - Jil||J,, 

where (■,■) is the inner product, and Z, Z\, and Z 2 are Lagrange multipliers. With 
appropriate fixed Z, Z 1; and Z 2 , an optimal C , D, L, R can be found by simply 
minimizing E(C, D, L, R, Z, Zi, Z 2 ) with respect to (C, D, L, R). Therefore, the trick 
is to determine Z, Zi, and Z 2 . In the augmented Lagrangian algorithm, we use 

' ( C (*+D > D (*+i) > j^+U, ^ fc+1 )) = arg min £(C, £>, X, Z« zf >, zf ), 

C ,D,L, Ft 

Z (h+i) = Z (k) + ( V ( L (k+i) R (k+i)^ _ F ^ 
z ( k +i) = z (k) + (c(fe+ i) _ VL (k+i)^ 

Zf +1) = Zf > + - ^ fe+1 )). 

As for the first minimization problem in (jHJ), we use an alternative minimization 
algorithm. As such, we minimize L sequentially in four steps, each with respect to 
one of the variables C, D, L, and R at a time, while keep the rest fixed. These 
four subproblems can be solved easily by either entry-wise soft-thresholding (for C and 
D) and linear equation solvers (for L and R), e.g. conjugate gradient method. The 
algorithm is summarized in Algorithm HJ where T is the soft-thresholding operator 
defined by [T a (A)]ij = sign([A]y) • max{|[A]y| — a, 0}. 

Algorithm 1 

1: Solve the first subproblem in (Q by alternative minimization iteratively as follows 

(i) Update C <- Ti /tXl {VL - zf%0- 

(ii) Update D <- T\j m {TR - Z { 2 k) /fi 2 ). 

(iii) Update L <- arg min^ %\\V(LR) -F+ Z«/^||| + f\\VL -C- zfV/^lll- 

(iv) Update R ±- axgimn R %\\V(LR) - F + Z^ /n\\ F + !f\\FR- D - Z { 2 k) /^ 2 \\ 2 F . 

2: If \\V(L^ k+r> R^ k+r> ) - F||| is small enough, then return. 
3: Z( fc+1 ) = Z« + (V(L^R^) - F). 
4: Zf' +1) = Z[ k) + (C( fc+1 ) - VL^). 
5: Zf +1) = Zf + (D( fc+1 ) - TR^). 
6: k = k + 1 and goto Step 1. 



One complexity of this problem comes from the non-convexity of the objective 
function in ([6]). Since there may exist local minima, we have to make sure our algorithm 
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get a desired solution. Yet, in our numerical experiments, this non-convexity issue is 
not found to be a problem for the following two reasons. First of all, our algorithm 
may find the global minimum of ((6]). As pointed out in Burer and Monteiro (2005), the 
augmented Lagrangian algorithm attains the global minimum for non-convex objectives 
in low-rank factorization form under suitable assumptions. Similar algorithms to ours 
has been used in, e.g., SDPLR (Semi-definite programming via low rank factorization) 
(Burer and Monteiro, 2003). Secondly, we have also chosen the initial guess to the 
iterative algorithm in the Algorithm [1] carefully as following. We first solve a convex 
minimization problem 

min U\VU - F\\ 2 F + X\\U\U, (10) 
u 2 

where is the nuclear norm, i.e., the summation of the singular values of U. The 

nuclear norm minimization is able to find a lowest-rank solution t/* of ([5]) within a 
precision (Cai et al., 2010). Then, the initial guess in the Algorithm [1] is chosen to be the 
best rank- if approximation to [/*. Let [/* = WSV T be a singular value decomposition. 
If the rank of E7» is less than if, we choose L<® = W£ 1/2 and = Y, l l 2 V T 

as the initial guesses. If the rank of [7* exceeds if, we choose = Wk^k and 
^k^k as the initial guesses, where Wk, Vk are the first if columns of W and 
V respectively, and is the if x if principle submatrix of 5]. By this way, L^R^ 
is the best rank- if approximation of L7*. 

As for the problem in ([7]), we still use the Algorithm [fl but we stop the iteration as 
soon as \[P(LR) — F\\ 2 F < a 2 . By this way, we can get a quite good approximate solution 
to (Q. This approach has been previously used and discussed in similar mathematical 
problems but in other contexts, e.g., (Osher et al., 2005; Cai et al., 20096; Cai 
et al., 2009a). 

2.3. Choice of K 

One practical issue in our algorithm is the selection of the parameter if to control the 
matrix rank. In practice, we achieve this goal by the following procedure. 

We first set if to be a large enough number, e.g. if = 20, and run our algorithm. 
Recall that the columns of L forms a basis to represent the images in U. Because of the 
sufficiently large if value in this trial run, the algorithm is forced to generate a set of 
if basis vectors containing those good ones for U, as well as the unnecessary bad ones. 
This is reflected by the fact that some columns of L in the solution make significant 
contributions to U, while the other columns contribute little. For those unnecessary 
columns in L, their presence introduces a small but observable amount of signals into the 
finally reconstructed images in U and compromises its quality. It is therefore desirable 
to eliminate them from the first place by setting the value of K to be the number of 
the significant columns. A second reconstruction is then performed with this properly 
chosen K value. 

It remains to determine which columns are significant and which ones are negligible. 
Since the contribution of the i-th column of L to U is L(:,i)R(i, :), any norm of this 
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project reflects the importance of this basis vector L(:, i). The larger norm is, the more 
significant L(:,i) is. In the experiments, we use the oo-norm to identify the significance, 
which is the maximum of row absolute sums of a matrix. The column L(:, i) is considered 
to be negligible, if \\L(:,i)R(i, :)\\oo is close to zero, and significant otherwise. 

2.4- Experiments 

We tested our cine-CBCT reconstruction algorithm on a digital NURBS-based cardiac- 
torso (NCAT) phantom (Segars et a!., 2008), which generates a patient body in thorax 
region with detailed anatomical features and a realistic motion pattern. The patient 
respiratory period is 4 seconds. The CBCT gantry rotates about the patient at a 
constant speed for a full rotation in 59 seconds, in which 360 X-ray projections are 
acquired. At each X-ray projection acquisition, we compute the NCAT phantom 
image at the specific time point with a resolution of 128 x 128, and the projection 
is then computed using a ray-tracing algorithm at the associated projection angle with 
a detector resolution of 256 bins. The patient breathing period and the gantry rotation 
period are chosen to be incommensurate deliberately, so that all the 360 patient images 
are distinct, although some of them visually look close. Some of the underlying true 
images are shown in Figure [TJ Under this setup, the size of U is 16384 x 360, and the 
size of F is 256 x 360. 

Two experiments are performed in this feasibility study, namely with and without 
noise signals in the measured projections F. For the case with noise, a realistic CBCT 
projection noise signal corresponding to 0.5 mAs/projection level is generated according 
to a noise model (Wang et al., 2008) and is added to the noise-free measurements F. In 
each case, we will first present the results for the selection of the parameter K and then 
the reconstruction results. Apart from visual inspection of the reconstructed cine-CBCT 
images, we have also quantitatively assessed the restored image quality using relative 
error as a metric defined by \\LR — U tr ue\\F l\\Utrue\\Fi where L, R are the outputs of 
the proposed method, and U trU e is the matrix consisting of ground truth NCAT images. 

3. Experimental Results 

3.1. Noiseless data 

For this case, the plot of ||X(:, i)R(i, :)||oo in a descending order as a function of the 
column index % is first depicted in Figure |2} with a trial run of K = 20. Clearly, 
the 8th column and beyond contribute little to the reconstructed images and hence 
should be removed. As such, K = 7 is selected and the corresponding results are 
shown in Figure [3] and Figure HJ Comparing the restored images in Figure |3] with the 
corresponding ground truth images in Figure [U it is found that our algorithm is able to 
capture the motions of the anatomy and restore the structures, even those small ones 
inside the lung. Meanwhile, observable artifacts inside the heart and at its boundary 
also exist. Quantitatively, the relative error of the restored cine-CBCT images is 3.97%. 
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Figure 1. The ground truth cine-CBCT images. From left to right, top to bottom: 
frame 40, 80, ...,360. 
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Figure 2. ||£(:, i)R(i, :)||oo in descending order. The result is produced with K = 20. 
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Figure 3. The restored cine-CBCT image from noiseless data. The relative error is 
3.97%. From left to right, top to bottom: frame 40, 80, ... , 360. 



To further look into the reconstruction results, we plot the columns of L in 
Figure 11(a), where each column is reshaped into a 128 x 128 image. These images 
form the basis to represent all the reconstructed cine-CBCT images. It is observed 
that the first one is similar to an image averaged over all the cine-CBCT images. Its 
presence provides an overall structure that is common to all the images in the cine- 
CBCT. Meanwhile, other basis images represent differences between images of U. We 
also plot the corresponding coefficients in R in Figure E(b), which attain a periodically 
variation pattern, indicating the patient respiratory motion. 

3.2. Noisy data 

For the case with noise in the projection data, the parameter K is selected as K = 4 
according to Figure El Figure [6(a) depicts the reconstructed images. The relative error 
of the restored cine-CBCT images is 6.81%. The columns of L and the rows of R are 
plotted in Figure [6(b) and (c), respectively. Again, we see that the columns of L are 
meaningful basis corresponding to the average image as well as the variations between 
them, and the rows of R oscillate periodically. 
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(b) The rows of R. 



Figure 4. L and R in the cine-CBCT image reconstruction from noiseless data. 



4. Conclusion and Discussions 



In this work, we have proposed an innovative algorithm to reconstruct cine-CBCT 
to provide a set of time-dependent images in the thorax region. In contrast to the 
currently used 4DCBCT, where only respiratory phase-resolved imaging is achieved, 
the cine-CBCT modality offers a much higher temporal resolution, as the instantaneous 
patient anatomy based on a CBCT projection is retrieved. Cine-CBCT reconstruction 
is apparently a challenging problem due to the very limited projection information 
associated to each image to be reconstructed. Yet, by effectively incorporating 
the underlying temporal constraints satisfied by those cine-CBCT images into the 
reconstruction process, we demonstrate the feasibility of this modality. In particular, 
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Figure 5. i)R(i, :)||oo in descending order. The result is produced with K = 20 

and noise at 0.5 mAs/projection. 

we assume a matrix factorization form U = LR for the matrixe U whose columns 
are the images to be reconstructed. The dimension of the matrices L constraints the 
rank of U and hence implicitly imposing a similarity condition among all the images in 
cine-CBCT. Moreover, near periodicity of the breathing pattern is reflected in R, which 
is achieved by minimizing the sparsity of its Fourier coefficients. Simulation studies on 
an NCAT phantom serve as a very preliminary test to demonstrate the feasibility of our 
approach. The relative error for the reconstructed images is 3.97% for the noise-free 
case and 6.81% for the case with noise. 

In addition to the apparent advantage of improved temporal resolution in the cine- 
CBCT over 4DCBCT, the reconstruction of cine-CBCT is not as demanding regarding 
the amount of projection data as 4DCBCT, at least based on our simulation studies. 
Because of the phase binning in 4DCBCT, a certain number of projections are required 
in each bin to yield a CBCT image for the bin with a satisfactory quality. Hence, 
it is usually challenging to get a 4DCBCT with only one-minute scan and in practice, 
protocols of multiple rotations or a slow rotation is utilized. In contrast, the cine-CBCT 
effectively incorporates the inter-image correlations into the reconstruction and a set of 
images of acceptable quality are obtained even with a single-rotation scan protocol, as 
demonstrated in our studies. 

Besides the modality of cine-CBCT, another potential application of our approach 
is the patient breathing signal extraction based on projection images. The respiratory 
signal is important for many clinical applications, such as projection sorting for 4DCBCT 
reconstruction. Currently, this signal is usually obtained by monitoring a surrogate 
known to move in synchronization with respiratory motion, e.g. patient surface or 
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(a) The restored cine-CBCT image with the relative error 6.81%. 
From left to right, top to bottom: frame 40, 80, . . . , 360. 




(b) The columns of L. Each column is reshaped into an image. 




100 200 300 100 200 300 100 200 300 100 200 300 

(c) The rows of R. 



Figure 6. Result of cine-CBCT image reconstruction with noise at 0.5 
mAs/projection. 
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diaphragm. However, these surrogates may not be readily available for a particular 
patient and there also may exist phase shift between the motions of the lung and that 
of the surrogate. Therefore, it is desired to obtain a robust breathing signal from 
the projections itself. As pointed in our algorithm, the matrix entries of the matrix 
R contains the breathing information, as evidenced by Figure [3](c) and Figure EJ^c), 
especially for those signals with a relatively large amplitude. Hence, our algorithm 
provides an effective method to extract the breathing signals from only projection 
information. 

Despite the success in this feasibility study, there are a few practical issues as for 
the clinical applications of cine-CBCT. First, only one slice of the patient anatomy is 
reconstructed in this proof-of-principle study. When cone-beam geometry is considered 
to reconstruct volumetric CBCT images, although the fundamental principles in our 
algorithm can be easily applied, there potentially come some problems due to the 
extremely large data size involved. The computational efficiency will also be reduced. 
This problem may be resolved or alleviated by using more powerful computational 
platforms, such as computer graphic processing unit (GPU) (Jia et al., 2010). Yet, the 
algorithm implementation and optimization are not straightforward for the practical 
considerations of the limited memory and the non-trivial parallel computation schemes 
in many steps. For instance, solving the least square problems in the Algorithm [1] 
requires carefully implementations of matrix multiplications and the associated matrix- 
vector operations. These issues require further investigations. Another issue that 
may potentially deteriorate the image quality is the irregularity of breathing pattern. 
The penalty term regarding the sparsity of the Fourier coefficients in favors those 
breathing patterns with a high level of periodicity. While this approach effectively 
imposes the periodicity aspect of the temporal correlation among cine-CBCT images, 
its performance, however, may be degraded in realistic patient cases with occasionally 
irregular breathing motions, such as cough. We plan to modify our current model 
to specifically take occasional irregularity into account and perform further tests in 
phantom data with realistic breathing motion and in real patient cases. 
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